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Abstract: We have characterised Water-based Liquid Scintillator (WbLS) using low energy pro¬ 
tons, UV-VIS absorbance, and fluorescence spectroscopy. We have also developed and validated a 
simulation model that describes the behaviour of WbLS in our detector configurations for proton 
beam energies of 210 MeV, 475 MeV, and 2 GeV and for two WbLS compositions. Our results have 
enabled us to estimate the light yield and ionisation quenching of WbLS, as well as to understand 
the influence of the wavelength shifting of Cherenkov light on our measurements. These results are 
relevant to the suitability of WbLS materials for next generation intensity frontier experiments. 
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1. Introduction 

Water-based liquid scintillator (WbLS) is a recently developed scintillating material that has been 
identified as a candidate detector medium for the next generation of intensity frontier particle 
physics experiments OH]. WbLS is a stable scintillating emulsion with a large fraction of aque¬ 
ous phase. 

With suitable chemical purification and material optimisation, WbLS may be produced with 
an optical attenuation of tens of meters in the photon energy range relevant to detection by bial¬ 
kali photomultiplier tubes (PMTs) (figure |T]). Such a long opfical attenuafion length would make 
Cherenkov imaging feasible in a large WbLS detector, which is not possible in large detectors 
filled with ordinary liquid scintillator. As a scintillating material, WbLS is sensitive to particle 
interactions below the Cherenkov thr'eshold. 

The aqueous phase of WbLS may also allow the loading of metallic ions that are not easily 
incorporated into the non-polar organic solvents that are typical of liquid scintillators. This may 
permit novel schemes to incorporate metal ions at high concentrations in the WbLS material. We 
are currently evaluating the possibility to incorporate candidate isotopes for neutrinoless double 
beta decay, neutron absorbing isotopes, and high Z materials for gamma attenuation. 
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Figure 1. The optical attenuation spectra for the 1% WbLS used in this study (solid black) and for a 
more recent formulation of 1% WbLS (dashed black). The 1% WbLS fluorescence spectrum (green) and 
photomultiplier tube (Hamamatsu R7723) quantum efficiency (red) are also indicated. The fluorescence 
spectrum was measured using an exciting wavelength of 290 nm. 


The optical and chemical properties of WbLS suggest that it may be an effective detection 
medium to be used in a particle detector capable of a broad range of measurements. Neutrinoless 
double beta decay, high energy neutrino beam measurements, proton decay, and low energy neu¬ 
trino physics have been proposed as plausible measurements using a large scale (30 - 100 kiloton) 
WbLS detector |0]. 

In order to rigorously evaluate the suitability of WbLS for such an application, a comprehen¬ 
sive understanding of the performance of the material must be developed. Basic material properties 
such as the light yield, the optical attenuation length, and the emission spectrum must be measured. 
A model of the response of the WbLS to ionising radiation and optical and ultraviolet photons must 
also be developed and validated. 

In this study, we have performed several experimental tests aimed at elucidating these prop¬ 
erties, and we have also developed a simulation model that successfully predicts the performance 
of the WbLS material in our measurement geometries. It is our intention to further develop and 
validate this simulation model using larger volume geometries to permit accurate study of the per¬ 
formance of a proposed large experiment that uses WbLS as its detection medium. 

2. Experiment 

2.1 WbLS Samples 

WbLS samples composed of 0.4% (WbLS-1) and 1% (WbLS-2) liquid scintillator by mass were 
prepared, as well as pure liquid scintillator (LS) samples. The pure LS was identical to that used 
by the Daya Bay experiment |0]. Table |I] details the compositions of these samples. 

2.2 Optical Characterization Measurements 

Ultraviolet-visible absorption spectroscopy of water, WbLS-1, and WbLS-2 were taken using a 
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Table 1. The chemical composition of the scintillators used in this study. The organic phase is given as 
a mass fraction of the final solution. The fluor concentrations are given relative to the organic solvent. 
Chemical names are abbreviated as follows; PC = Pseudocumene, PPO = 2,5- diphenyloxazole, bis-MSB = 
1,4- Bis (2-methylstyryl) benzene, LAB = Linear Alkyl Benzene, LS = Liquid Scintillator. 


Sample Name 

Organic Phase 

Primary Fluor 

Secondary Fluor 

WbLS-1 

PC, 0.4% 

PPO, 0.4 g/L 

bis-MSB, 3 mg/L 

WbLS-2 

PC, 0.99% 

PPO, 1.36 g/L 

bis-MSB, 7.48 mg/L 

LS 

LAB, 100% 

PPO, 2 g/L 

bis-MSB, 15 mg/L 
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Figure 2. Excitation-emission spectrum for WbLS-2. 


Shimandzu UV-1800 spectrophotometer. The optical attenuation spectrum of WbLS-2 is shown in 
figure H], along with that of a more recent WbLS formulation that also consists of 1% scintillator, 
the fluorescence emission spectrum, and the photomultiplier quantum efficiency spectrum for the 
PMTs used in this study. The optical attenuation of the more recent formulation of a LAB-based 
WbLS is markedly less than WbLS-2, after the purification of all starting materials and the addition 
of a surface-scattering reduction agent. 

Fluorescence spectroscopy was carried out using a PTI fluoresence spectrometer. The excitation- 
emission spectrum of WbLS-2 is shown in figure which exhibits an excitation wavelength de¬ 
pendence of the emission spectrum shape. 

2.3 Proton Beam Measurements 

Two beamline instruments were developed for this study. Detectors A and B (figure 

Detector A was fabricated using black ABS polymer to supress light reflections from the walls. 
The 130 mm x 115 mm x 63.5 mm detector volume was filled with water during initial irradiations, 
then WbLS-2. Two PMTs (Hamamatsu R7723, 51 mm diameter) observed the liquid volume; 
one placed downstream from the beam so as to allow detection of Cherenkov light, and one placed 
upstream from the beam so as to suppress detection of Cherenkov light. To allow for any variability 
in the photomultiplier response, all measurements were repeated with the PMT locations relative 
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Figure 3. A schematic drawing of the instruments used in the proton beam measurements, Detector A and 
Detector B. See the text for details. 


to the beam direction swapped. The PMT response was found to be similar for each tube (figure 
We report the averaged values of the PMT response hereafter. Two 2x2x0.5 cm plastic scintillator 
hodoscopes (HI and H2) were used for triggering and to define the beam location. 

Detector B consists of two right cylindrical 150 mm x 150 mm tubs. Tub 1 was fabricated using 
6.35 mm thick white teflon and Tub 2 was fabricated using 6.35 mm thick aluminium coated with 
black perfluoroalkoxy alkane paint on its inner surface. The tubs were exposed to the proton beam 
for consecutive fillings of water, WbLS-1, WbLS-2, and LS. Two 2x2x0.5 cm plastic scintillator 
hodoscopes were used to define the beam. 

A 400 nm LED was placed just above each scintillating volume in Detector A to permit an in- 
situ single photoelectron PMT calibration. The low energy water irradiations were used for single 
photoelectron PMT calibrations in Detector B. 

All samples were investigated using proton beams with incident energies of 475 MeV and 2 
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Figure 4. The measured photoelectron distributions for identical measurements using the two different pho¬ 
tomultipliers (represented as the red and blue traces) in Detector A. The solid traces represent measurements 
taken downstream of the 2 GeV proton beam on a water target, and the dashed traces represent measure¬ 
ments taken upstream of the 475 MeV proton beam on a WbLS-2 target. The photomultiplier responses are 
similar for small and large numbers of photoelectrons. 


Table 2. A summary of the proton beam exposures. 


Instrument 

Sample 

Incident Proton Energy 

Detector A 

Water 

2 GeV and 475 MeV 

Detector A 

WbLS-2 

2 GeV and 475 MeV 

Detector B 

Water 

210 MeV, 475 MeV, and 2 GeV 

Detector B 

WbLS-1 

210 MeV, 475 MeV, and 2 GeV 

Detector B 

WbLS-2 

210 MeV, 475 MeV, and 2 GeV 

Detector B 

LS 

210 MeV, 475 MeV, and 2 GeV 


GeV at the NASA Space Radation Laboratory at Brookhaven National Laboratory. A 210 MeV 
incident proton beam was additionally used to investigate all samples in Detector B. The beam in¬ 
tensity at all energies was made low enough that the probability of having more than a single proton 
incident per accelerator bunch was small. An analysis of our results where the two-proton events 
are separable from the 1 proton events indicates that the fraction of two-proton events varied as a 
function of beam energy and was less than 10% in all irradiations. We have estimated the influence 
of the two-proton events by resampling the simulated single proton distributions presented below. 
Double proton events have no effect on the results with a large (>10) mean number of detected 
photoelectrons. The effect on the remaining data is estimated to be <10% and is included in the 
systematic uncertainty discussed in section ^ 

The data readout in all measurements was achieved using a CAEN V1729A 14 bit waveform 
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Figure 5. A schematic of the waveform data processing algorithm used to determine the pulse information 
for analysis. The dashed red line is the baseline determined by the algorithm. See the text for details. 

digitizer sampling at 1 gigasamples per second. The raw waveforms were stored for offline analysis. 
Data acquisition was triggered by coincident hodoscope events during the accelerator beam gate. 

3. Analysis 

3.1 Signal Processing Algorithm 

Signal processing of the acquired waveforms was achieved using a custom algorithm which was 
developed for background subtraction, timing, and pulse area determination. The algorithm is 
based upon that used in the DarkSide experiment 0] and its operation is illustrated in Figure ^ 

A moving average of the waveform data is taken as the baseline. A data point is consid¬ 
ered to be part of the baseline if the difference between that data point value and the value of the 
averaged baseline as it was trig_pre_samples data points ago is less than the algorithm’s 
threshold parameter, max_amplitude. If a trigger event occurs (trig_start), the algo¬ 
rithm waits until the waveform falls below the threshold level, then waits for an additional pile-up 
pulse for untrigger_length samples. If a pile-up event occurs, trig_stop is unset and 
the piled-up event is considered to be part of the same trigger. Otherwise, the algorithm linearly 
interpolates between the baseline values at the start and end of the pulse and continues to average 
the baseline and search for a new trigger. The pulse charge was measured as the summed pulse 
area between trig_pre_samples and trig_post_samples, less the baseline. We used 
trig_pre_samples = 3 ns and trig_post_samples = 10 ns. 

The use of a moving average of the baseline attenuates high frequency components of the 
signal pedestal noise; allowing a max_amplitude value of approximately 0.3 photoelectrons. 
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3.2 Simulation Model 


We have used Geant4 vlO.O [01 to simulate detectors A and B, the energy deposit in the hodoscopes 
and scintillator, and the optical photon measurements by the PMTs. 

For modeling the optical photon processes, the optical properties of the detector materials 
were required. Where possible we used measured values. Our measurements of the excitation and 
emission spectra and optical attenuation coefficient spectrum of the scintillators used in this study 
were incorporated into the model. The wavelength-shifting absorption coefficient was taken to be 
equal to the difference between the WbLS optical attenuation coefficient and the optical attenua¬ 
tion coefficient of water. The scintillation emission spectrum was assumed to be identical to that 
produced when excited by 290 nm photons. The refractive index of water and optical absorbance 
of water were taken from Segelstein [0. The refractive index and re-emission probability of liquid 
scintillator were taken from the Daya Bay simulation model [0]. The refractive index of WbLS 
samples was calculated as a linear combination of water and pure liquid scintillator, according to 
the fraction of scintillating solvent present in the WbLS. The spectral dependence of the WbLS 
re-emission probability was assumed to be identical to that of pure liquid scintillator, although the 
amplitude of the re-emission probability was left as a free parameter whilst being constrained to 
a maximum value of 1. The 2 GeV measurements in Detector A were employed for this optimi¬ 
sation. We assumed that the re-emission probabilities of WbLS-1 and WbLS-2 are identical. The 
optical absorbance and refractive index data for the UV-transparent acrylic windows were taken 
from Band et al. [H]. The refractive index of the borosilicate glass PMT windows was calculated 
using the empirical Sellmeier equation parameters provided by Schott [0. The quantum efficiency 
spectrum was estimated using data provided by Hamamatsu. 

Several optical parameters were either unknown or poorly known, so that their values were as¬ 
sumed or left as free parameters. The black vessels were assumed to be perfectly unreflective, and 
the white vessel was assumed to be a wavelength-independent diffuse reflector with the reflectance 
allowed to vary. The water measurements were used to calibrate other model parameters to repro¬ 
duce the measured results. The free parameters were the PMT photocathode radius for Detector A 
and the beam height relative to the centre of the detector vessel for Detector B. The beam height in 
the Detector B was estimated to be located at the centre of the tub with an uncertainty of ± 1 cm. 
The optimal beam height was found to be 0.69 ± 0.05 cm above the centre of the vessel. The op¬ 
timal photocathode radius of 22 ± 1 mm is in fair agreement with the manufacturer’s specification 
of >23 mm. 

3.2.1 Wavelength Shifting Model 

For ordinary liquid scintillator with a scintillation yield of 10"* photons per MeV, the number of 
scintillation photons produced per centimeter for a minimum ionising particle experiencing a stop¬ 
ping power of 2 MeV/cm is 2 x 10"*. This vastly outnumbers the yield of Cherenkov photons, 
which is Ri 600 photons per centimeter in water integrated over the 200-500 nm spectral range. 
As more than half of these photons fall in the 200 to 300 nm spectral range, the Cherenkov light 
accounts < 1.5% of the detectable light from a minimum ionising particle in ordinary liquid scintil¬ 
lator. WbLS has a much lower scintillation light yield than ordinary liquid scintillator, so that the 
number of Cherenkov photons may in fact be greater than the number of scintillation photons. This 
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Figure 6. The predicted optical photon spectrum in Detector A due to the scintillation (dashed traces) and 
wavelength-shifting (solid traces) processes, as measured by the upstream PMT for incident 2 GeV protons. 
The original Geant4 wavelength-shifting model results are represented using red traces, and the results from 
the modified wavelength-shifting model developed in this work are represented using blue traces. 


Cherenkov light may be directly detected or absorbed by the WbLS. Wavelength-shifting (WLS) of 
the absorbed light may occur with a probability described by the WbLS’s re-emission probability. 
Understanding the Cherenkov and scintillation light generation and WLS is therefore an important 
aspect of evaluating the suitability of WbLS for large detector geometries. 

This sensitivity to WLS of Cherenkov photons has motivated us to extend the Geant4 optical 
WLS model. Our modifications to the WLS model have allowed both the re-emission probabil¬ 
ity and the emission spectrum to depend upon the wavelength of the absorbed photon. Figure ^ 
compares the model predictions of the optical photon spectrum incident on the upstream PMT in 
Detector A due to a beam of 2 GeV protons, when it is filled wifh WbLS-2. Allhough a large 
proporlion of wavelengfh-shiffed lighf is predicfed by bolh models, Ihe original Geanl4 model pre- 
dicls 21% less WLS lighf. There is also some specfral disforlion evidenf in fhe modified WLS 
model specfrum af around 450 nm due fo fhe influence of fhe changes in fhe fluorescence emission 
speclrum lhaf occur mainly for exciting wavelengfhs befween abouf 350 nm and 400 nm (figure 0). 
Our modified WLS model predicfs a mean number of defecled phofons lhaf is 17% grealer fhan fhe 
Geanf4 WLS model. 

3.2.2 Scintillation Parameter Optimisation 

The scintillation light yields and ionisation quenching factors of all samples were taken as free 
model parameters. The light yields were estimated by performing a optimisation of the sim¬ 
ulated photoelecton distribution to the measurement, using data from Detector B when irradiated 
by 475 MeV protons. The ionisation quenching factor optimisation was performed in a similar 
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Figure 7. Measured (red traces) and simulated (blue traces) photoelectron distributions for the scintillator 
beam tests in Detector B, Tub 2. 


manner, using Detector B’s 210 MeV proton irradiation. 

4. Results 

4.1 Detector B 

The simulated and measured photoelectron (PE) distributions are given in figures 0 and Cuts 
were applied to the hodoscope pulse amplitude and timing, as well as the relative PMT-hodoscope 
timing to select single proton events from the incident beam. The cuts effectively removed the 
pedestal noise and did not bias the photoelectron distribution. The measurement of 475 MeV 
and 210 MeV protons on a water target are not shown as few events were observed during this 
measurement, and those events that were registered were of very low amplitude. We attribute the 
events observed in these measurements to the detection of Cherenkov emission from secondary 
electrons. The mean number of observed photoelectrons for all measurements and simulations 
are given in table H. In addition to the statistical uncertainty arising from the analysis of the 
simulated data, we estimate that there is an additional systematic uncertainty component of ss 10% 
with the main contributions due to uncertainties in the simulation model’s input data (such as the 
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Figure 8. Measured (red trace) and simulated (blue trace) photoelectron distribution for a 2 GeV proton 
beam incident upon water in Detector B, Tub 2. The range over which the minimisation to the measured 
data was performed is shown. 

Table 3. The most probable energy deposit and mean number of photoelectrons (N/*^) in Detector B, Tub 
2. The most probable energy deposit is taken from the Monte Carlo simulation. We estimate a systematic 
uncertainty of « 10% for the simulated values. Where uncertainties are given as 0.0, the actual uncertainty 
is smaller than the smallest reported significant figure. 


Sample 

Proton Beam Energy 

Most Probable 
Energy Deposit 

NpE, Measured 

NpE, Simulated 

Water 

210 MeV 

106.4 ± 0.3 MeV 

1.1 ±0.0 

1.5 ±0.2 

Water 

475 MeV 

42.0 ± 0.2 MeV 

1.7 ±0.0 

1.3 ±0.1 

Water 

2000 MeV 

26.6 ± 0.3 MeV 

2.8 ±0.1 

2.5 ± 0.3 

WbLS-1 

210 MeV 

113.9 ±0.3 MeV 

4.5 ±0.1 

4.4 ± 0.4 

WbLS-1 

475 MeV 

42.2 ±0.1 MeV 

3.6 ±0.1 

3.4 ± 0.3 

WbLS-1 

2000 MeV 

27.4 ± 0.2 MeV 

17.0 ±0.2 

18.6 ± 1.9 

WbLS-2 

210 MeV 

113.6 ±0.3 MeV 

27.6 ± 0.3 

27.3 ± 2.7 

WbLS-2 

475 MeV 

42.1 ±0.2 MeV 

16.7 ±0.1 

15.9 ± 1.6 

WbLS-2 

2000 MeV 

27.5 ± 0.3 MeV 

30.8 ± 0.3 

29.2 ± 2.9 

LS 

210 MeV 

96.6 ± 1.4 MeV 

2588 ±21 

2622 ± 262 

LS 

475 MeV 

32.5 ± 0.6 MeV 

nil ± 10 

1105 ± no 

LS 

2000 MeV 

20.9 ± 0.4 MeV 

872 ± 15 

933 ± 93 


refractive index, PMT quantum efficiency, and re-emission probability), the simplified treatment 
of optical photon scattering, and the effect of two-proton events (Section ^3| ). This systematic 
component dominates the uncertainty of all simulated results. 
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For Tub 2, the simulated results are generally in agreement with the measurements. The 
optimal values of light yield and ionisation quenching parameter are given in table ^ We have used 
Birks’ semi-empirical ionisation quenching model [HQ] , and the quenching parameter refers to the 
material parameter of that model. The light yield of the LS sample is consistent with the typically 
measured value of the light yield of a liquid scintillator of ss 10^ photons/MeV. The quenching 
parameter for the LS sample falls in the lower range of the typically measured values of 0.07-0.2 
mm/MeV [O o, n Q. The light yield values for the WbLS samples do not linearly scale 
with the fraction of scintillator, likely due to differences in fluor concentration and solvent type. 
Both WbLS samples are more susceptible to ionisation quenching than the LS sample. Indeed, the 
quenching parameter for both WbLS-1 and WbLS-2 is larger than any previously reported value 
that could be found by the authors for a liquid scintillator. The cause of these extraordinary values 
will be investigated in future studies. 

The 475 MeV and 210 MeV measurements in WbLS-1 resulted in small numbers of photo¬ 
electrons with broadened distributions relative to a Poisson distribution. It is also apparent that 
some low amplitude events with « 1 PE contributed to the distribution, as evidence of these events 
are seen in some of the other measurements. The distribution broadening can be modeled by con¬ 
volving the simulated distributions with a Gaussian as outlined in [I13ll . However the broadened 
distribution obtained using this technique were unable to reproduce the shape of the measured pho¬ 
toelectron distributions for all proton energies. Due to the poor agreement of distribution shapes, 
optimisation was quite sensitive to the selection or exclusion of the 1 PE bin in the analysis. 
The sensitivity of the light yield and quenching parameter to this effect is included in the statistical 
uncertainties for the WbES-1 results reported in table 0. 

The Tub 1 results are not presented, as it was not possible to obtain agreement with the sim¬ 
ulation model that was consistent with the measured distributions across all beam energies for all 
materials, for any value of the tub reflectance. While the cause of this inconsistency is unknown, 
we have identified two possible causes that arise from simplifications in our simulation model. 
Eirstly, there may be some wavelength-dependence to the tub reflectance. If this is the case, the ef¬ 
fective reflectance for Cherenkov light is different from the effective reflectance of the scintillation 
and WLS emission, which may explain the discrepancies seen with different particle energies and 
samples, that have different proportions of Cherenkov and scintillationAVLS light. Another source 
of error that may contribute to the discrepancy is the model’s treatment of optical attenuation. In 
the model, the optical attenuation length - the sum of optical absorption and scattering - is treated 
as equivalent to the optical absorption length. Treating scattering as equivalent to absorption is not 
generally a problem in an unreflective detector with low PMT coverage - a photon that is scattered 
is likely to be absorbed at the wall. Therefore we do not expect this mechanism significantly in¬ 
fluence the results in Tub 2. However, in a highly reflective detector such as Tub 1 the scattered 
photons are more probable to be detected, so that the conflation of scattering with absorption is 
a source of systematic error. This issue will need to be addressed prior to the application of the 
model to large detector geometries with greater PMT coverage. 

4.2 Detector A 

The simulated and measured photoelectron (PE) distributions for Detector A are given in figure ^ 
Quality cuts to the data were applied in the same way as Detector B. The measurement of 475 MeV 
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Table 4. The optimal values of light yield and ionisation quenching parameter for each scintillator. 
Material Light Yield (photons/MeV) Quenching Parameter (mm/MeV) 


WbLS-1 19.9 ±1.1 (stat.) ± 2.0 (sys.) 0.70 ± 0.12 (stat.) ± 0.07 (sys.) 

WbLS-2 108.9 ± 0.8 (stat.) ± 10.9 (sys.) 0.44 ± 0.01 (stat.) ± 0.04 (sys.) 

LS 9156 ±42 (stat.) ±916 (sys.) 0.07 ± 0.01 (stat.) ± 0.01 (sys.) 


Table 5. The mean number of measured and simulated photoelecttons (N/)^) for Detector A, for photomul¬ 
tipliers located downstream and upstream of the proton beam. We estimate a systematic uncertainty of « 
10% for the simulated values. Where uncertainties are given as 0.0, the actual uncertainty is smaller than the 
smallest reported significant figure. 


Sample 

Incident Proton Energy 

Photomultiplier 

NpE, Measured 

NpE, Simulated 

Water 

475 MeV 

Downstream 

1.3 ±0.0 

1.6 ±0.2 

Water 

475 MeV 

Upstream 

1.2 ±0.1 

1.2 ±0.1 

Water 

2000 MeV 

Downstream 

33.0 ±0.2 

32.4 ± 3.2 

Water 

2000 MeV 

Upstream 

1.2 ±0.0 

1.4 ±0.1 

WbLS-2 

475 MeV 

Downstream 

4.7 ± 0.0 

4.6 ± 0.5 

WbLS-2 

475 MeV 

Upstream 

4.6 ± 0.0 

4.5 ± 0.5 

WbLS-2 

2000 MeV 

Downstream 

21.5 ±0.3 

20.4 ± 2.0 

WbLS-2 

2000 MeV 

Upstream 

7.7 ± 0.2 

7.3 ±0.7 


protons in water are not shown. The mean number of observed photoelectrons for all measurements 
and simulations are given in table ^ 

The simulated results used the optimised values of WbLS-2’s light yield and quenching from 
the Detector B analysis. Whilst fairly good agreement is obtained using these values, the shape of 
the simulated photoelectron distribution is in general broader than that of the measured distribution. 
The cause for this behaviour is not fully understood. We have also performed a optimisation 
of the simulated distribution to the measured one by adjusting the simulated light yield using the 
475 MeV proton irradiation data. These results suggest an optimal light yield of 109.0 ±1.0 and 
110.2 ±1.0 photons per MeV using the downstream and upstream PMTs, respectively. The un¬ 
certainties relate only the statistical uncertainty arising from the x^ analysis of the simulated data, 
and we estimate an additional 10% systematic uncertainty. This good agreement with the Detector 
B analysis suggests that the model can consistently describe different measurement geometries. 

For the 2 GeV proton measurements, many more photons were observed in the downstream 
PMT than the upstream PMT due to the directional Cherenkov radiation. However, the number 
of photoelectrons measured in the downstream PMT at 2 GeV for WbLS was substantially less 
than the number of photoelectrons measured in water at the same energy, which suggests that some 
proportion of the detectable Cherenkov light is absorbed or scattered by the WbLS. Due to the black 
detector walls, the upstream detector is largely insensitive to the direct Cherenkov light (as shown 
in figure 0^b)), and therefore this signal can be used to estimate the number of photoelectrons that 
are measured due to the isotropic light emission processes - scintillation and WLS. The absorption 
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Figure 9. Measured (red traces) and simulated (blue traces) photoelectron distributions for the beam tests 
using Detector A. 

probability of the detectable Cherenkov light in the WbLS is given to first order by: 

P I 2 GeV)-Nup{WbLS ,2 GeV) 

Ndown{Water,2 GeV) 

where Ndown and Nup are the mean numbers of photoelectrons in the downstream and upstream 
photomultipliers for a given measurement, respectively. Our measurements indicate that Pats ~ 
58%; that is, 58% of Cherenkov photons that were measured in the water sample were instead 
absorbed in the WbLS sample measurement. 

The large proportion of Cherenkov photons absorbed or scattered by the i^(lOcm) path length 
in WbLS suggests that the optical attenuation length of the WbLS used in this measurement is 
too short to permit Cherenkov imaging in a very large detector. The optical attenuation length of 
more recent WbLS formulations has been improved by about 2 orders of magnitude (figure |I]). To 
firsf order, fhe fraction of absorbed Cherenkov phofons simply scales wifh fhe opfical affenualion 
lengfh, which suggesfs fhaf fhe ss 58% loss of Cherenkov phofons would occur over a pafh lengfh 
of (^(lOm) using fhe more recenf WbLS formulation; which is a lengfh scale suifable for a large 
defector. Experimenfal measuremenfs in a larger defecfor geomefry are needed fo more rigorously 
assess fhis model prediction. 

Our measuremenfs can also give an esfimafe of fhe relative proporfions of measured WLS and 
scinfillafion phofons. The mean number of phofons defecfed in fhe upsfream PMT of Defecfor A 
for 2 GeV profons incidenf on WbLS-2 may be wriffen as fhe sum of fhe defecfed scinfillafion and 
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WLS light components: 

Kp{WbLS,2 GeV) = Nf,f\WbLS,2 GeV) +Nup^^^^''‘’'’\wbLS,2 GeV) (4.2) 

where N^p"^ {WbLS, 2 GeV ) and {WbLS, 2 GeV ) are the mean number of detected scintil¬ 

lation photons and WLS Cherenkov photons, respectively. Any wavelength shifting of scintillation 
light is disregarded in this approximation. 

The mean number of photons detected in the upstream PMT of Detector A for 475 MeV 
protons incident on WbLS-2 may be written as a similar expression, though without any WLS 
Cherenkov light component, as 475 MeV protons are at the Cherenkov production threshold: 

NupiWbLS,425 MeV) = Nl‘j;"\WbLS,415 MeV) (4.3) 

The 2 GeV and 475 MeV protons experience stopping powers in water of 2.0 and 2.8 MeV/cm, 
respectively. It is possible estimate the difference in scintillation light yield at the different energies 
using the Birks quenching correction: 

dE 

m) = ( 4 . 4 ) 

where E is the proton energy and Ub is the quenching parameter. Note that equation assumes a 
constant stopping power for the particles interacting in the tub. The corrected light yields can then 
be used to relate the number of scintillation photons measured using the 475 MeV proton beam 
with the 2 GeV proton beam: 

Nup{WbLS,2 GeV) = -^f^^fupiWbLSAlS MeV) +Nup^^^^""\wbLS,2 GeV) (4.5) 

Equations permit the calculation of the relative contributions of the scintillation 

and Cherenkov processes to the measured signal in the upstream PMT. Our results indicate that 
in the upstream PMT for incident 2 GeV protons; for every detected photon that was produced by 
scintillation 1.27 ±0.05 photons are detected that are produced by WLS Cherenkov light. This 
value is consistent with the simulation prediction of 1.28. The relatively high proportion of re¬ 
emitted Cherenkov light suggests that the detection sensitivity of a WbLS-based detector will be 
improved signihcantly beyond what would be expected by the scintillation light yield alone for 
particles that exceed the Cherenkov threshold in the medium. However, the WLS Cherenkov light 
is less localised to the position of the energy deposit in the detector and modihes the linearity of the 
energy response of the scintillator. These complicating factors will need detailed study in a larger 
measurement geometry in order to evaluate their effect upon track reconstruction and calorimetry. 

5. Summary and Outlook 

The work outlined in this study has made a hrst measurement of the light production and prop¬ 
agation characteristics of water-based liquid scintillator. Two different WbLS concentrations and 
pure liquid scintillator were studied. The simulation model developed in this work appears to be 
fairly robust as it was able to reproduce the measured photoelectron distributions in kilogram-scale 
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detectors for three different proton energies that each probed different light emission properties 
of the WbLS. Higher quality WbLS re-emission probability measurements and better treatment of 
scattering will be used to improve the simulation model. Measurements and model validation in 
larger test geometries are also required before our model can be used with confidence to predict the 
performance of a kiloton-scale detector. We have undertaken to develop a 1000 liter scale WbLS 
detector for this purpose. 
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